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Our conclusions J,J , about the value of p in the time 
evolution to oc i^ of the magnetization, are questioned 
in Ref. Q- In order to dispel the suspicion that our 
Monte Carlo (MC) results may be size dependent, we 
show in Fig. 1 results for various system sizes, both in 
SC and FCC lattices, that (1) exhibit that p does depend 
on lattice structure, as reported in Ref. [ij, and (2) show 
no trace of any size dependence. We find size effects only 
in smaller systems. 

The variation of p with lattice structure should not be 
perplexing. It is predicted by our theory [2J, of which 
the main result, i.e., Eq. (11) of Ref. [l|, is known. Its 
predictions, given in Ref. [l|, are also shown in Fig. 1. 
Furthermore jfor vanishing tunneling window widths, Eq. 
(11) of Ref. Q gives 



p ho 



(1) 



for all times well after F"^, where a is approximately the 
dipole field rms value |3], ho is Stt^/S^/^ times a nearest 
neighbor dipole field, and F"^, as well as all other no- 
tation in this Reply, are as defined in Ref. m. Values 
for a are given in Ref. |3 for various lattices, but, for 
completeness' sake, we give here cr/ho = 0.756, 0.398, 
0.417, and 0.66 for SC, BCC, FCC, and Fcg lattices, re- 
spectively. Values for p that follow from Eq. (1) agree 
rather well with the MC results we have reported and 
with experiments on Feg |^. 

So why are the numerical results of the Ref. [21 so 
different from ours? Ours yj apply to (1) annealed 
pI systems, to which a magnetic field H is applied at 
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FIG. 1: m, scaled with ew£aH/a^ , versus Vt for SC and FCC 
lattices with N = L x L x L and A x L x L x L/2 spins, 
respectively. All points stand for averages over at least 10^ /N 
MC runs, with £„ = 0.05, H = 1, and £a = 0.43 and 0.50, 
for SC and FCC lattices, respectively. Full lines stand for 
theoretical results. 



ksT < O.IU/S, if (2) i/ < cr [such as iJ < 4 mT and 
cr ss 30 mT (i.e., gfisS^oH < 0.05 K and gpsS^oO' ~ 0.4 
K in obvious notation) in Ref. ^]. Neither of these two 
conditions is met in Ref |2|, where initial spin configura- 
tions are random and H > a. 

Annealing is essential in the magnetization process 
studied in Ref. [l| . The very nature of the process de- 
pends on it. In annealed systems, spin-up and spin-down 
populations that are able to tunnel, i.e., on sites where 
dipolar fields approximately cancel iJ, are unequal, and 
this drives the magnetization process. In unannealed 
systems both populations are, on the average, initially 
equal, and much slower thermal equilibration processes 
then drive the magnetization evolution j^. Furthermore, 
annealing is hard to avoid in experiments such as in Ref. 
y. An Feg crystal held for as little as 1 second within 
the range 20 > T > 2 K before quenching to much lower 
temperatures, will qualify as annealed [5j, and its magne- 
tization will be at least two orders of magnitude larger, 
for up to minute (roughly 4F~^) after H <^ a is applied, 
than it would have been had the initial spin configuration 
been somehow completely randomized initially. 

Furthermore, the results we have obtained for the mag- 
netization of annealed systems also apply to the relax- 
ation of the magnetization in zero field, after cooling in 
a weak field Q. Then, m(0) - m{t) ex t^ while 1 < Ft 
and m{t) > 0, and p is also given by Eq. (1). 

Incidently, in unannealed systems (as in Ref. [1]), p is 
also nonuniversal, varying with H , increasing monotoni- 
cally as H decreases, up to a value slightly larger than 1 
for H <^ a. 
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